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Low  Order  Adaptive  Optical  System 

Abstract 

The  report  represents  the  evaluation  results  of  the  performance  of  the 
bimorph  piezoelectric  adaptive  mirror.  The  first  part  of  results  concerns 
spatial  modes  of  the  mirror,  their  linearity,  hysteresis  and  their  temperature 
dependence.  The  next  part  of  results  concerns  the  dynamic  behaviour  of  the 
mirror. 

The  last  part  of  the  report  represents  results  of  the  performance  of  the 
simple  adaptive  optical  system  for  imaging  through  laboratory  generated 
turbulence.  In  conclusion  we  presented  the  general  evaluations  of  the  adaptive 
optical  system  with  bimorph  mirror  for  lm  class  telescope. 
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1.  Bimorph  piezoelectric  deformable  mirror. 


1.1.  The  progress  in  adaptive  optics  (AO)  application  is  mainly 
associated  with  increasing  the  number  of  controlled  subapertures  and 
feedback  control  circuits.  As  a  result  of  such  development  there  appeared 
systems  with  more  than  200  subapertures  and  the  plan  was  reported  to  built 
490  subapertures  AO  telescope  [1].  There  are  two  major  problems  on  this  way. 
First,  it  leads  to  increasingly  high  cost  of  AO  systems.  Second  -  the  higher 
the  order  of  correction  required  (this  broadly  equals  the  number  of  the  mirror 
actuators),  the  more  light  you  need  to  get  accurate  measurements  of  the 
distorted  wavefront. 

One  way  to  overcome  these  problems  is  to  use  a  rather  small  number 
of  active  subapertures  for  on-line  real  time  pre-processing  of  the  observed 
image  and  to  optimise  parameters  of  AO  system  to  achieve  effectiveness 
required  for  post-processing. 

In  experiments  on  AO  phase  correction,  performed  in  MSU  since 
1979,  we  used  bimorph  deformable  mirrors  (BDM)  with  comparatively 
small  number  (from  8  to  17)  of  actuators  or  spatial  deformation  modes  as  key 
elements.  A  number  of  papers  were  published  on  specific  problems  concerning 
such  DM  applications: 

-  parameters  optimisation  and  design  [2]; 

-  measurements  of  operating  parameters  [3]; 

-  computations  of  effectiveness  in  different  adaptive  systems  [4]. 

A  number  of  laboratory  experiments  were  performed  on  BDM 
applications.  Results  obtained  revealed  sufficiently  high  effectiveness  of 
bimorph  mirrors  to  compensate  low  order  aberrations  such  as  thermally 
induced  distortions  in  laser  systems  [5,6].  The  last  reports  by  other  groups 
show  the  wide  interest  for  BDM  astronomical  application  as,  for  example  - 
COME  ON  PLUS  and  CFH  projects  [6]. 
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1.2.  The  concept  of  BDM  has  been  suggested  in  [7]  and  at  first  realised 
by  Lipson  and  Steinhaus  [8],  The  device  (see  fig  1.1)  consists  of  a  glass  plate 
firmly  glued  to  a  plate  actuator  disc,  made  from  piezoelectric  ceramic.  In  our 
BDM  the  actuator  disc  on  its  turn  consists  of  two  piezoceramic  discs  (lead 


zirconium  titanate),  soldered  together  and  polarised  normally  to  their  surfaces. 

Maoism  _ 


The  diameter  of  the 
glass  plate  is  39  mm  and 
thickness  4  mm.  The  active 
surface  of  the  plate  has 
reflecting  aluminium  layer 
with  protecting  SiC>2 


cover.  The  thickness  of 
each  piezoceramic  disc  is 
0.35  mm. 

The  interface 


Fig.  1.1  Bimorph  adaptive  mirror 


between  the  two  discs 


contains  continuos  conducting  electrode  -  we  shall  call  it  "ground  electrode". 
The  second  continuos  conducting  electrode  between  the  actuator  and  the  glass 
plate  is  used  to  control  the  curvature  of  the  whole  mirror  and  we  shall  call  it 
"general  focus  electrode".  The  outer  surface  of  the  actuator  has  8  separate 
ring-segment  electrodes  to  control  mirror’s  astigmatism. 


The  mirror  profile  dependence  via  control  voltages  applied  to  electrodes 
can  be  described  by  Poisson  equation  [8,9]. 
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1.3.  In  the  experiments  presented  we  used  two  mirrors  (fig.  1.2).  One 
mirror  was  fixed  to  the  mount  by  a  glass  tube  attached  to  the  actuator  disc. 
We'll  call  it  BDM1  -  "fixed  centre  mirror".  The  second  mirror  was  fixed  to 
the  mount 
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Fig. 1.2.  Mirror  attachment  variants: 
fixed  centre  (left),  quasi-free  mirror 


by  a  soft  plastic  bandage.  We'll  call  it  BDM2  -  "quasi  free  mirror".  The 
BDM1  mount  has  no  tip-tilt  actuators  so  we  should  use  separate  plate  mirror 
for  tip-tilt  control.  Actually  it's  not  a  disadvantage,  because  any  modem 
telescope  has  a  separate  precise  tracking  system.  Moreover,  the  fixed  centre 
mirror  can  be  built  in  the  AO  set-up  more  precisely,  than  quasi  free  mirror. 


The  BDM2  hold  mount  can  be  attached  to  the  tip-tilt  controlled  mount. 
This  mount  has  2  PZT  stick  actuators  with  2.5  pm  stroke  at  300  V  drive  that 
allows  us  to  get  about  5  mrad  tip-tilt  control  range. 

Due  to  the  large  capacitance  of  actuators  and  the  relatively  high  weight 
of  the  attached  BDM2  holder  mount  the  dynamic  parameters  of  the 
configuration  are  rather  poor  -  it  has  less  than  100  Hz  frequency  bandwidth. 
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The  tip-tilt  controlled  mount  can  be  used  for  angular  correction 
separately  by  attachment  of  a  small  (up  to  25  mm)  plate  mirror.  This 
configuration  allows  to  obtain  1.2  KHz  bandwidth  that  is  quite  enough  for 
astronomical  turbulence  compensating  systems. 
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2.  Spatial  parameters  of  the  mirror 


Practically  all  results  presented  in  this  section  concerns  the  BDM1  - 
fixed  centre  mirror.  Spatial  parameters  of  the  BDM2  have  been  completely 
measured  previously  in  MSU.  Shortly  we  can  summarise  that  the  second  order 
aberrations  may  be  fitted  with  RMS  error  about  6%  and  max  p-p  aberration 
value  is  4  pm  at  300  V.  The  first  resonance  frequency  was  about  1-1.2  KHz 
and  active  bandwidth  is  about  500-600  Hz. 

2.1.  Interferometer.  Spatial  parameters  of  mirrors  were  measured  on 
commercial  digital  phase  shifting  Moller-Wedel  interferometer  with  Phase 
Shift  Technology  software.  The  main  operating  parameters  of  the 
interferometer: 

Phase  measurement  resolution  -  A./256  (^=0.63pm); 

Phase  measurement  accuracy  better  than  V50  p-p; 

Repeatability  better  than  V100; 

Data  acquisition  time  less  than  1  second; 

Fitting  of  phase  data  to  36  Zemike  polynomials. 

Data  obtained  on  the  interferometer  can  be  stored  in  three  file  formats: 
INT  -  for  intensity,  MAP  -  for  phase  profile,  and  ZRN  -  for  Zemike 
representation.  Application  software  to  convert  this  data  into  standard  BMP 
and  ASCII-DAT  file  formats  was  written  using  Borland  Turbo  Pascal  7.0.  The 
most  part  of  results  is  presented  in  terms  of  Seidels  aberrations. 

To  control  BDM  aberrations  we  used  the  handle  control  unit,  that 
consists  of  6  potentiometers  and  2  DC  voltage  converters  -  from  +12V  to 


+300V  and  from  +12V  to  -300V  DC.  Output  voltages  of  this  unit  can  varies 
from  -150  V  to  +150  V. 


2.2.  Low  order  aberrations  fitting.  Hysteresis. 

2.2.1.  Plane  fitting.  The  interferogram  and  spatial  profile  of  the  mirror 
with  zero  control  voltages,  measured  at  18°C,  are  shown  at  fig.2.1.  The  main 
initial 
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Tilt  Focus  Astigm  Coma 
Sphercl 

0.12  -1.35  0.24  0.05 

-0.10 

TOTAL 

Peak:  0.68  Valley:  -0.82 
P-V:  1.50 

RMS:  0.394  (Seidel:  99.86%) 
Strehl:  0.000 


Fig.2.1.  The  profile  of  the  mirror  at  zero  control  voltages 

aberrations  are  general  focus  -  1.35  pm  (p-v),  astigmatism  -  0.24  pm,  and 
spherical  0.1  pm,  the  last  one  is  mostly  due  to  the  centre  fixing  rod  influence. 
Average  RMS  value  of  aberration  was  0.394  pm,  and  the  most  part  of  it  - 
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99.86%,  was  represented  by  Seidel  aberrations:  defocus,  astigmatism,  pure 
coma  and  spherical. 

We  connected  the  opposite  segmented  electrodes  of  the  mirror  so  we 
applied  one  control  voltage  to  el&e5  electrodes  pair,  the  next  control  voltage 
to  e2&e6  pair  and  so  on.  By  handle  control  of  these  5  voltages  (one  for  general 
focus  electrode)  we  succeeded  to  obtain  plain  fitting  with  0.169  pm  p-p  and 
0.033  pm  RMS  aberration  -  fig.2.2.  Control  voltages  for  this  case  are 
represented  in  Table  2.1. _ _ 


u 


_ X/R _ 

p-p  ABERRATIONS  in  pM 

Tilt  Focus  Astigm  Coma 
Sphercl 

0.07  0.067  0.022  0.089 

-0.165 

TOTAL 

Peak:  0.098  Valley:  -0.071 
P-V:  0.169 

RMS:  0.033  (Seidels:  91.8%) 
Strehl:  0.899 


Fig.2.2.  Plane  fitting  by  5-channels  handle  control 


Table  2.1.  Control  voltages  for  plane  Fitting. 


el=e5 

e2=e6 

e3=e7 

e4=e8 

e9 

5.6 

-9.2 

-2.3 

16.4 

-42 
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The  maximum  voltage  -42V  was  applied  to  electrode  e9  (general  focus 
electrode)  to  compensate  initial  defocus  of  the  mirror.  For  the  computer 
controlled  closed  loop  adaptive  system  with  a  wave  front  sensor  we  may  hope 
to  obtain  the  fitting  of  plane  with  less  RMS  error. 

2.2.2.  General  focus  control.  The  dependance  of  the  Seidel  aberrations 
via  the  closed  cycle  of  the  e9  control  voltage  variation  are  represented  at 
fig.2.3  (general  focus)  and  2.4.(astigmatisms,  coma  and  spherical  aberrations). 
The  curve  at  fig.2.2  shows  high  controllability  of  the  mirror  by  general  focus 
electrode:  the  amplitude  of  general  focus  variation  6.28  jam  at  245  V  control 
voltage  variation.  These  figures  gave  us  the  next  value  of  the  control  gain 
coefficient  kp: 

kp  =  AF/AV  =  0.0256  jiM/V. 
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These  results  reveals  the  typical  influence  of  the  PZT  ceramic  hysteresis, 
for  general  focus  aberration  the  absolute  value  of  the  width  AFj2of  the 
hysteresis  curve  was  0.56pM  which  corresponded  to  the  relative  value  8p  = 
8.9%.  In  closed  control  circuite  the  influence  of  hysteresis  can  be  significantly 
decreased  by  additional  feedback  loop  [12].  Control  gain  coefficients  for  Seidel 
aberrations  are  shown  in  table  2.2. 


Fig.2.3.  Seidel  aberrations  for  general  focus  (e9)  control. 


Table  2,2.  Control  gain  coefficients  for  e9  general  focus  control. 


Aberration 

coma 

I1H55I 

control  gain,  xlO^  pM/V 

-6.1 

-3.6 

5.7 

The  similar  curves  (see  fig. 2. 4,  2.5)  were  obtained  for  general  focus 
control  by  applying  the  same  voltages  to  all  segmented  electrodes  -  el..e8.  As 
it  could  be  expected  the  control  gain  coefficient  kp  was  almost  2  times  less  in 

compare  with  the  previous  case  (control  by  electrode  e9): 

=  AF/AV  =  0.01 12  pM/V. 


n 


Control  gain  coefficients  for  Seidel  aberrations  for  segmented  electrodes 
control  are  shown  in  table  2.3. 


Aberration 

coma 

control  gain,  xl()4  pM/V 

-3.4 

-0.57 

16.7 

2.2.3.  Spherical  aberration.  The  significant  difference  in  the  control  gain 
coefficients  of  the  spherical  aberration  for  e9  and  el-e8  cases  gives  us 
possibility  to  control  this  aberration  by  applying  a  certain  voltage  difference 
between  e9  electrode  and  the  ring  of  segmented  electrodes.  The  results  for 
spherical  aberration  control  are  shown  at  fig. 2. 6  -  spherical  aberration  and 
fig.2.7  -  astigmatism,  coma  and  focus.  At  fig.2.6  horizontal  axis  represents 
segmented  electrodes  voltage  and  corresponding  e9  voltage  values  are  shown 
near  experimental  points.  Control  gain  coefficient  for  spherical  aberration  had 
the  value  11 0'3  mM/V. 
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0.4, 


Negligible  values  of  general  focus  aberration  (fig. 2.7)  proves  that  we 
imposed  correct  differences  between  e9  and  el-e8  voltages.  The  significant 
variations  of  astigmatism  in  this  case  can  be  explained  by  a  certain 
misalignment  between  the  centre  of  the  mirror  plate  and  the  centre  of  fixation 
rode  as  well  as  by  imperfection  of  electrodes  geometry. 

2.2.4.  Astigmatism.  To  control  astigmatism  we  applied  a  certain  positive 
voltage  to  a  pair  of  opposite  segmented  electrodes  and  the  same  negative 
voltage  to  another  pair  of  electrodes,  as  it  shown  at  fig.2.8.  To  compensate  the 
initial  defocus  we  applied  20.2V  to  the  e9  electrode.  Axes  of  maximum 
deformation  for  this  case  almost  coincided  the  X-Y  axes  and  in  terms  of 
Seidels  we  can  call  this  aberration  2°-astigmatism.  The  control  gain  coefficient 
for  this  case  had  the  value  £ast=2.2-10'2  pm/V. 
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Fig.2.8.  Astigmatism  fitting  results. 


The  similar  results  for  reversed  signs  of  control  voltages  gave  us  the  value  of 
£ast=1.95  10'2  (j.m/V.  To  obtain  45°-astigmatism  we  applied  +139V  to 

electrodes  e4,e8  and  -139V  to  electrodes  e2,e6.  The  control  gain  coefficient 
£as j.  for  this  case  had  the  value  1.88  10'2|im/V  and  2.23  10‘2  pm/V  for 

reversed  signes  of  applied  voltages. 
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2.2.5.  Coma.  Originally  this  mirror  with  8  segmented  electrodes  wasn't 
designed  to  correct  the  coma  aberration.  As  it  was  difficult  to  fit  the  coma  by 
handle  control  of  6  voltages  available  the  next  procedure  was  used.  For  every 
segmented  electrode  we  measured  the  response  profile  and  their  Zemike 
expansions  (36  terms)  were  used  to  calculate  the  accuracy  of  the  lowest 
aberrations  fitting  by  LSE  method.  The  results  of  these  calculations  are 
represented  in  table  2.4  (fitting  errors)  and  2.5  (control  voltages). 


Table  2.4.  Calculated  Errors  of  the  lowest  Zernike  Fitting 


Aberration 

P-to-V,  pM 

RMS,  pM 

RMS/P-V,  % 

Plane 

- 

0.0315 

Focus 

1.27 

0.0076 

1.27 

0.037 

2.9 

1.27 

0.038 

3.1 

Pure  Coma  1 

1.27 

0.06 

4.8  ! 

Pure  Coma  2 

1.27 

0.2 

15 

Spherical 

0.949 

0.28 

30 

Table  2.5. 

Calcul 

ated  Control 

Volta 

?es  to 

fit  the 

lowes 

Zemike,  V 

Aberration 

el 

e2 

e3 

e4 

e5 

e6 

e7 

e8 

e9 

Plane 

RH | 

-20 

-6.8 

16.4 

16.7 

7. 

-7.1 

8.2 

8 

Focus 

5.9 

6 

6.1 

5 

5.5 

6 

10.7 

3.9 

-42.2 

QSSSS3DHH 

31.4 

5.3 

-28.8 

0.7 

33.9 

5 

-24.8 

3.8 

3.8 

M4 1 W  lii— 

0.5 

31.8 

■El 

-32.9 

mm 

26.5 

3.8 

-30.2 

Pure  Coma  1 

10.2 

-86.8 

-132 

-89.3 

-34.2 

54.4 

80 

64.1 

-19.7  | 

Pure  Coma  2 

30.4 

24.9 

9.2 

-16.4 

-28.8 

-15.3 

-2.7 

22.1 

■u 

5.4 

9.3 

15.7 

9.8 

1 

-5.1 

-7.3 

■W 

In  accordance  with  the  results  presented  in  tab.2.4  while  one  type  of  the 
coma  -  "pure  coma  1"  can  be  fitted  with  RMS  error  less  than  5%,  another  one 
-  pure  coma  2  has  significantly  higher  fitting  error  -  15%  in  term  of  RMS. 
The  fitting  error  value  obtained  for  spherical  aberration  is  about  30%.  At  the 
same  time  fitting  errors  for  the  second  order  aberrations,  for  which  this 
mirror  was  designed,  are  about  3%  for  astigmatism  and  0.6%  for  defocus. 

2.2.6.  Estimations  for  turbulence  correction  application.  For  general 
estimation  of  the  mirror  effectiveness  in  astronomical  AO  system  we  can  use 
results  obtained  by  R.J.Noll  [9].  He  calculated  weights  of  Zemike  aberrations 


16 


in  wavefront  distortions  caused  by  Kolmogorov's  type  turbulence.  These  low 
aberrations  weights  in  the  total  phase  variance  on  a  circular  aperture  D 

are  represented  in  table  2.6.  In  the  table  2.7  the  figures  of  the  vision  quality  in 
terms  of  the  Strehl  coefficient  S  are  represented  for  a  number  of  D/ro  values 
(ro  -  wavefront  coherence  length). 

Table  2.6.  Low  order  aberrations  percentage  in  the  total  phase  variance  c£on 


circular  aperture  (exl.pision) 


Tilt 

Focus 

Astigmatisms 

Comas 

Spherical 

residual 

87% 

2.23% 

2.23%  x  2 

1.2%  x  2 

0.23% 

3.66% 

Table  2.7  Vision  quality  for  different  atmospheric  conditions 


D/r0 

<r;=1.03-(D/r„r 

S»e^in 

s, 

compensated 

T  +  F  + A 

s, 

compensated 

T  +  F  +  A  + 

C  +  Sph 

1 

1 .0299 

0.598 

0.968 

0.981 

2 

i 

3.27 

0.194 

0.902 

0.942 

5 

15.057 

. -  . 

5  10-4 

0.623 

0.753 

10 

47.81 

4  10-H 

0.222 

0.417 

This  figures  show  that  compensation  of  the  second  order  aberrations 
can  significantly  improve  the  vision  quality  for  small  (1  -  1.5  m)  telescopes  at 
average  values  of  ro  about  10-20  cm.  For  the  same  values  of  D/ro  we  can 

calculate  the  desired  deformations  of  adaptive  mirror  surface  to  correct  these 
second  order  aberrations  -  see  table  2.8.  In  this  calculations  we  used  the 
relation  Dx=3o<!>  that  means  that  for  normal  distribution  of  aberrations  in 
95%  of  time 

wavefront  distortion  p-p  value  will  be  in  the  range  (-3ad>..  +3ad>). 

Table  2.8.  Desired  profile  deformatons  in  micrometers  (for  95%  probability) 


D!ra 

1 

2 

5 

10 

Ax  =  -  3??. -;(A  =  2 ptM) 

In 

0.97 

1.73 

3.71 

6.6 

Ax  TILT 

0.9 

1.61 

3.45 

6.14 

Ax  FOCUS,  ASTIGMATISM 

0.15 

0.26 

0.56 

0.99 

Ax  PURE  COMA 

0.11 

0.19 

0.41 

0.73 

Ax  SPHERICAL 

0.05 

0.08 

0.18 

0.32 

The  greatest  desired  value  is  about  1  pm  for  the  second  order  aberration 
deformation.  The  same  is  for  coma.  These  results  enable  us  to  hope  that  the 
mirror  investigated  can  be  used  successfully  to  compensate  atmospheric 
turbulence. 


2.3.  Temperature  stability.  The  temperature  stability  of  the  mirror 
parameters  is  of  great  importance  for  the  real  applications.  In  the  room 
where  interferometer  is  installed  (it  is  small  enough  )  we  use  the  simple  heat 
gun  for  increasing  the  temperature  in  the  range  of  a  few  °C  degrees.  For 
temperature  measurements  digital  thermocouple  type  thermometer  Thandar 
TH302  was  used  (precision  0.1  °C).  The  thermocouple  was  attached  to  the 
mirror  mount. 

At  fig. 2.9  the  interferogram  and  topography  of  the  mirror  are  shown 
for  temperature  values  18°C  and  22.6°C.  This  figures  shows  significant 
temperature-defocus  dependence  due  to  bimorph  structure  "glass  plate  - 
piezoceramic  actuator  disk"  as  well  as  to  the  difference  in  temperature 
expansion  coefficients  of  the  glass  and  piezoceramics  which  were  used  to 
fabricate  the  mirror. 
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Fig.  2.9.  The  interferogram  and  topography  of  the  mirror  for  2  temperature 


values 


The  temperature  dependencies  of  aberrations  in  terms  of  p-v  values  are 
represented  at  fig. 2. 10.  To  match  the  experimental  results  the  linear  LSE 
approximation  was  used.  The  most  significant  value  for  temperature 
dependence  coefficient  was  measured  for  defocus  -  0.38  jim/°C,  while 
coefficients  for  others  low  order  aberrations  are  1 0-30  times  less. 
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Astigmatism:  Coma:  0.001 1*T  Spherical:  0.0025-T  - 

0.003  T  +  0.04  +0.008  0.047 

Fig.2.10.  Temperature  dependencies  of  the  lowest  aberrations. 

As  we  can  see  from  results  presented  in  previous  sections,  the  range  of 
aberration  control  for  this  mirror  is  high  enough  to  correct  turbulence  induced 
distortions  and  to  compensate  its  own  thermal  distortions  at  10-15°C  variations 
of  temperature. 
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3.  Temporal  parameters. 


The  resonance  frequencies  of  the  mirror  was  measure  with  the  help  of 
Kenwood  AG-204  oscillator,  ISO-Tech  ISR-620  oscilloscope  and  Black  Star 
Apollo- 10  frequency  meter.  During  the  measurements  we  connected  driving 
sine  voltage  to  X  input  of  the  oscilloscope  as  well  as  to  frequency  meter  while 
the  signal  measured  was  gone  to  the  Y  input  of  the  oscilloscope.  This  enabled 
us  to  use  both  phase  and  amlitude  measurements  of  resonances  that  was  more 
precise  than  pure  amplitude  measurement. 


Fig. 3.1. 


The  frequency  response  curve  shown  at  fig.3.1  was  obtained  by  applying 
AC  voltage  to  e9  electrode  (general  focus),  output  signal  was  measured  on  el 
electrode.  Similar  curves  for  the  case  when  driving  voltage  was  applied  to  the 
el  electrode  and  response  was  measured  at  one  of  another  segmented 
electrodes  are  shown  on  fig.3.2. 
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The  lowest  resonant  frequency  of  the  mirror  is  about  6  KHz,  while  small 
peak  at  2.5  KHz  -  (see  fig  3.2c)  may  be  caused  by  the  resonance  of  the  mirror 
attachment  plate. 


frequency,  KHz  frequency,  KHz  frequency,  KHz 


a  b  c 

Fig. 3.2 

We  can  conclude  that  the  mirror  has  operating  frequency  range  up  to 
about  1  KHz,  what  is  sufficiently  enough  for  astronomical  applications. 


4.  Low  order  adaptive  system  performance 


4.1.  Experimental  set-up.  In  order  to  illustrate  the  operation  of  the 
BDM  we  constructed  a  simple  adaptive  system  for  imaging  trough  laboratory 
generated  turbulence.  Experimental  set-up  of  the  low-order  adaptive  system  is 
shown  at  fig.4.1.  The  beam  of  He-Ne  laser  with  10  mW  output  was  expanded 
from  1mm  to  10mm  width  by  lenses  LI  and  L2  and  directed  into  the  glass 
water  cell,  that  was  used  to  produce  turbulent  distortions  of  the  wavefront.  The 
turbulence  was  produced  by  convection  movement  of  water  between  the 
bottom  hot  metal  plate  and  the  cold  one  placed  at  the  top  of  the  cell.  After 
plane  mirror  Ml  one  part  of  the  beam  went  through  the  lenses  L5,L6  at  the 
BDM  surface.  The  beam  width  at  the  mirror  surface  was  35mm.  After 
reflection  this  beam  was  directed  in  the  "adaptive"  measuring  arm  where  lenses 
L3(F=200mm),  25  pm  pin-hole  D3  and  photodiod  PD]  were  used  to  measure 
Strehl  criterion  of  the  corrected  beam. 


Fig.4.1.  Adaptive  system  experimental  set-up. 


Another  part  of  the  laser  beam  was  gone  by  the  beam  splitter  directly 
into  measuring  arm,  containing  the  lens  L4(F=200mm),  25  |xm  pin-hole  D4 
and  photodiod  PD2,  that  were  used  to  observe  beam  quality  -  axis  intensity 
without  phase  correction.  In  order  to  get  the  measured  criteria  -  PD]  and  PD2 
signals,  the  PCL812PG  card  was  used.  The  ADC  circuit  on  this  board  has  8 
analogue  inputs,  12  bits  range  with  programmable  input  range  0.3125,  0.625, 
1.25,  2.5  and  5  V  (amplitude).  The  maximum  rate  of  ADC  conversion  is  25 
KHz.  It  also  allowed  us  to  check  the  program  and  hardware  delay  in  the  circuit 
photodiod  -  ADC  -  computer  -  mirror  with  onboard  clock  oscillator  with 
accuracy  0.5  msec. 

For  mirror  control  we  used  a  DAC  unit  with  1 1  output  signals  and  DC 
Amplifiers  with  150V  amplitude  output.  The  computer  used  to  control  the 
experimental  set-up  was  Dell  PC  486  type,  66  MHz,  8  Mb  RAM.  The  software 
for  our  experiment  was  written  on  Borland  Turbo  Pascal  7.0  and  gave  us 
possibility  to  get  and  to  observe  the  PD!  and  PD2  signals,  to  control  the  mirror 
and  store  the  data  on  hard  disk.  We  used  the  simplest  "hill-climbing"  [11] 
algorithm  to  calculate  our  control  signals  for  maximisation  of  the  measured 
criterion. 

4.2.  Turbulence  in  the  water  cell.  Experimental  data  for  two-point 
coherence  function  of  spherical  in  the  water  cell  were  obtained  by 
A.Bogaturov  and  V.Myakinin  from  Institute  of  Physics  of  atmosphere, 
Moscow  in  the  following  way. 

A  point  source  of  spherical  wave  (focused  laser  beam)  was  installed  at 
distance  1  m  from  the  turbulent  cell  and  at  distance  2  m  from  the  receiving 
lens.  (To  obtain  coherence  function  for  a  plane  wave  it  is  necessary  to  divide 
the  distance  by  2).  Average  images  of  this  source  were  captured  in  the 
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conjugated  plane  in  the  absence  of  turbulence  and  in  presence  of 
turbulence.  The  window  of  CCD-device  was  wide  open  and,  consequently, 
these  images  were  the  integrals  of  intensity  distribution  over  the  length  of 
one  pixel  (2.5  |i.m)  rather  than  intensity  distributions.  This  approach 
allowed  us  to  obtain  coherence  function  using  one-dimensional  fast  Fourier 
transformation  (FFT).  The  plot  of  this  coherence  function  is  shown  on  fig.4.2. 


0.0  0.2  0.4  0.6  0.8  1.0 


Ac  mm 

Fig.4.2.Two-point  coherence  function  of  spherical  wave; 

At  -  temperature  difference  between  cold  and  hot  plates. 

Assuming  that  the  two-point  coherence  function  may  be  presented  in  the 
following  form:  C(r)  =  expMr/ro)0*),  and  using  least-squares  method,  we  may 

conclude  from  these  data  that  a  approximately  equals  to  1 .9  for  all  three  cases 
of  temperature  difference. 
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In  our  experiments  with  adaptive  system  the  temperature  difference 
between  cold  and  hot  plates  was  30-35°C  and  tq  was  approximately  0.2  -  0.5 

pm. 


Fig.4.3. 

The  Criterion  dependencies  on  a 
separate  aberration  variation; 


The  Order  of  Control  voltage  variation: 
0V  ...  -150V  ...  +150V  ...  0V 


total  512  points  for  each  curve 
40  msec  point-to-point  delay 


4.3.  Controllability.  In  order  to  estimate  the  controllability  of  the 
system  at  the  first  stage  we  measured  the  influence  of  the  BDM  profile 
variations  on  the  beam  quality  criterion  -  PDj  signal  (with  no  turbulence). 
These  curves  for  X,Y,  defocus  and  astigmatisms  control  voltage  variations  are 
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shown  at  fig  4.3.  On  These  curves  reveals  the  strong  influence  of  tip-tilt 
actuators  hysteresis,  that  is  much  less  for  defocus  and  astigmatism  cases.  We 
can  consider  these  curves  as  2-D  cross  sections  of  the  beam  quality  criterion  in 
5-D  space  with  X,Y,Df,Astl  and  Ast2  coordinates. 

4.4.  Correction  of  static  aberrations.  At  the  second  stage  of 
experiments  we  checked  the  closed  loop  system  performance  for  correction  of 
static  aberrations  -  with  no  turbulence. 


Fig.4.4.  The  correction  of  the  static  aberreations. 


First  the  optical  set-up  was  completely  adjusted  and  the  value  of  our  PD] 
signal  was  measured.  Then  some  static  aberrations  were  imposed  by  small 
longitudinal  and  transverse  displacements  of  the  measuring  diaphragm  D3.  So 
the  diaphragm  was  slightly  shifted  across  the  focal  spot,  by  this  way  we 
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simulated  tip  and  tilt,  and  along  the  beam  axis,  that  gave  us  defocus 
aberration.  Actually  due  to  a  certain  misalignments  in  the  optics  we  obtained 
not  only  tip-tilts  and  defocus  but  some  others  aberrations  as  well.  After  that  the 
system  started  to  correct  the  imposed  aberrations.  As  we  didn't  need  high 
speed  performance  we  set  20  msec  step-time. 

The  typical  time  dependence  of  the  PDj  signal  during  adaptation  process 
is  shown  on  fig.4.4.  As  we  might  expect  the  maximum  value  of  the  criterion 
achieved  by  adaptation  varied  as  we  controlled  the  different  number  of 
aberrations.  So  the  final  PD]  value  for  only  tip-tilt  correction  was  108,  for  tip- 
tilt  and  defocus  -  211,  and  for  tip-tilt,  defocus  and  astigmatisms  correction  - 
228.  Initial  PDj  value  (before  aberrations  was  imposed)  was  225  and  it  was  48 
after  aberrations  were  produced  (just  before  adaptation  process  started). 

As  we  can  see  from  these  data  it  took  the  system  from  about  50  steps 
(for  tip-tilt  case)  to  about  1 50  steps  (for  5  aberrations  control)  to  achieve  the 
maximum  of  the  criterion.  The  small  variations  of  the  PD]  signal  after  the 
system  had  reached  the  maximum  were  due  to  permanent  "checking  up" 
variations  of  the  control  signals  in  accordance  with  the  "hill-climbing"  method 
[11]. 


4.5.  Correction  of  turbulence.  The  temporal  and  spatial  parameters 
of  the  turbulence  produced  in  the  water  cell  strictly  corresponded  each  other  - 
when  the  coherence  length  tq  was  large  the  typical  value  of  time  scale 

variations  was  large  as  well,  and  when  we  had  the  swift  turbulence  the 
coherence  length  was  small.  The  defraction  length  for  turbulence  induced 
phase  distortions  ro 2/X  for  ro=0.5mm,  X=0.63pm  had  the  value  about  40  cm. 

It  means  that  phase  distortions  in  the  mirror  plane  (as  well  as  in  the  measuring 
plane)  had  been  transformed  into  intensity  distortions  and  couldn’t  have  been 
completely  corrected. 


28 


No  Adaptation  Experiment  #133  j 

Tip-Tilt  Adaptation  Experiment  #:  134 

Cycle  Time:  2.6  msec 

Cycle  Time  :  2.7  msec 

<PDj>  =  89; 

<PD2>  = 

<PD]>=  107  <PD2>:  = 

93 

22. 

Strehl:  0.39 

Correlation:  0.352 

0.41 

Strehl:  0.47  0.098 

No  Adaptation  Experiment#:  135  Tip-Tilt  Adaptation  Experiment#:  136 


Cyrclc  Time  :  1 .6  msec  Cyrclc  Time:  1 .6  msec 

<PD|>  =  29  <PD2>  =  16  <PDj>:  76  <PD2>:= 

Slrehl:  0.13  0.065  13 

Correlation:  0.271 _ Strehl:  0.33 _ 0-053 

Fig.4.5  Tip-tilt  correction  of  laboratory  induced  turbulence; 

PD!  -  signal  after  Water  Cell  +  Adaptive  mirror  (solid  curves) 

PD2  -  signal  after  Water  Cell  (dashed  curves) 

_  PD]=PD2  =  228  without  turbulence _ 


The  typical  step  time  delay  in  the  system  caused  by  the  software  and 
hardware  delays  was  about  0.6  msec,  but  to  prevent  the  resonance  exitation  of 
the  tip-tilt  actuators  we  imposed  additional  programmable  delay  in  the  range  1- 
2  msec. 
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The  performance  of  the  system  in  the  presence  of  turbulence  is 
illustrated  on  fig.4.5  (tip-tilt  correction  only),  4.6  (tip-tilt+defocus  correction) 
and  4.7  (all  second-order  aberrations).  Left  columns  of  these  charts  show  the 
time  dependence  of  the  PD]  and  PD2  signals  (in  ADC  code  units)  in  opened 
loop  system,  and  the  right  columns  -  in  closed  loop  system. 


No  Adaptation  #147 
Cycle  Time:  1.6  msec 

<PD]>  =  44  <PD2> 

108 

Strehl :  0.19 
0.48 

Correlation:  0.454 


Tip-Tilt  +  Defocus  +  Astigmatisms 
#148 

Cycle  Time :  1.7 

<PD1>=173  <PD2>  = 

130 

Strehl  :  0.76 
0.57 


No  Adaptation  #149 

Cyrclc  Time  :  1.6  msec 

<PDj>  =  12  <PD2>  = 

77 

Strehl  :  0.005 

0.34 

Correlation:  0.194 

Tip-Tilt  +  Defocus  +  Astigmatisms 

#150 

CyrcleTime:  1.6 

<PD!>=161  <PD2>  = 

92 

Strehl:  0.71 

0.40 

Fig.4.7  Tip-tilt,  Defocus  +  Astigmatisms  correction  of  laboratory  induced 

turbulence. 

PD]=PD2  =  228  without  turbulence 

Corresponding  values  of  the  Shtrehl  factor  can  be  calculated  by  the  next 
formula: 

S 12  =  <PDj  2>/PDo, 
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where  we  denoted  sample  mean  value  by  brackets  <...>,  and  PDq  is  the 
maximum  value  of  criteria  without  turbulence.  In  the  described  experiments 
we  had  PDq=228. 


As  we  can  see  in  closed  loop  cases  the  average  values  of  the  "corrected" 
beam  criterion  -  PDj  were  higher  than  that  of  the  "uncorrected"  beam 
criterion  -  PD2.  The  comparative  values  of  the  Sj  and  S2  as  well  as  their 
standard  deviations  G!  and  a2  are  shown  in  the  table  4.1. 


Table  4.1. 


Corrected 

Aberrations 

uncontrolled  channel 

Adaptive  channel 

sjs. 

oJo2 

&  =<m  > 

a2 

5i  =<  PD,  > 

c. 

134  -  2ms 

tilts 

0.10 

0.10 

0.47 

0.20 

4.70 

2 

136 

tilts 

0.06 

0.06 

0.33 

0.17 

5.85 

2.8 

138 

tilts 

0.06 

0.07 

0.32 

0.15 

5.26 

2.1 

140 

tilts 

0.12 

0.10 

0.39 

0.15 

3.35 

1.5 

142 

tilts+defocus 

0.29 

0.18 

0.55 

0.26 

1.89 

1.44 

144 

tilts+defocus 

0.40 

0.15 

0.71 

0.08 

1.77 

0.53 

146 

tilts+defocus 

0.49 

0.13 

0.73 

0.05 

1.49 

0.38 

148 

0.57 

0.10 

0.76 

0.05 

1.33 

0.50 

150 

0.55 

0.11 

0.76 

0.05 

1.39 

0.45 

The  average  increase  in  the  Strehl  coefficient  obtained  during  the 
adaptive  system  performance  in  our  experiments  was  in  the  range  1.5-6.  We 
should  note  that  all  experiments  mentioned  in  table  4.1  could  be  subdivided 
roughly  in  two  groups.  In  the  first  group  -  mainly  tip-tilt  correction  -  we  had 
rather  great  Strehl  values  increase:  from  3.5  to  6,  while  in  the  second  group  - 
with  all  second  order  aberrations  corrected  (besides  tip-tilt)  -  we  had  only  1.4 
-  1 .5  increase.  At  the  same  time  in  the  first  group  we  had  2-3  increase  of  the 
standard  deviations,  while  in  the  second  group  the  criterion  variance  failed  in 
comparison  to  that  of  the  uncontrolled  channels. 


The  increase  of  the  criterion  standard  deviation  values  in  the  first  group 
can  be  explained  by  the  sharp  dependence  of  the  criterion  measured  on  tip-tilt 
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changes  -  see  fig.4.3.  It  means  that  we  had  exceedingly  large  control  gain  for 
tip-tilt  channels  that  caused  large  variances  in  criterion  even  at  small  control 
signals. 

Comparatively  small  increase  in  Strehl  values  in  the  second  group  in  its 
turn  was  due  to  the  increase  in  the  cycle  time  of  the  system  -  as  we  could  see 
from  part  4.4  the  typical  cycle  time  for  complete  aberration  control  was 
approximately  3  times  greater  than  that  of  the  only  tip-tilt  control  cases.  The 
decrease  in  the  o  values  in  this  group  showed  that  the  weight  of  the  second 
order  aberrations  in  the  beam  received  was  high  enough.  It  enabled  us  to  hope 
that  if  we  had  fast  measurements  of  the  beam  phase  profile  we  could  have 
much  better  results  on  adaptive  correction  with  our  mirror. 

The  next  step  of  this  system  development  could  be  the  additional 
complete  analogue  or  analogue-digital  circuit  for  tip-tilt  correction  with  rather 
simple  quadrant  detector  for  measurements  of  the  lens’s  L3  focal  spot 
movement.  The  final  stage  would  be  the  adaptive  system  with  Shark-Hartmann 
wave  front  sensor  with  relatively  small  number  of  light  sensitive  subapertures 
-  about  16x16. 
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5.  Conclusions 


The  main  spatial  parameters  of  the  mirror  investigated  are  shown  in  the 
table  5.1.  Here  we  can  see  the  control  coefficient  k,  the  temperature 
dependence  coefficient  j u,  relative  rms  of  aberration  fitting  and  p-p  aberration 
deformation.  The  plane  fitting  rms  is  about  0.033p.M. 


Table  5.1.  Main  parameters  of  the  mirror. 


Aberration 

k,  hM/V'102 
exprmnt 
calc 

at  200  V 

RMS/P- 

V 

% 

Focus 

2.56 

5.1 

0.6 

0.38 

Astigmatism 

2.06 

4.1 

3 

0.003 

Pure  Coma 

0.51 

1 

5  -  15 

0.0011 

Spherical 

0.1 

0.2 

30 

0.0025 

The  typical  width  of  the  mirror’s  actuator  hysteresis  curve  is  about  10% 
while  operating  frequency  is  almost  lKHz. 

We  had  built-in  the  mirror  in  the  simple  hill-climbing  adaptive  system 
for  compensation  of  the  laboratory  induced  turbulence.  The  performance  of 
this  system  showed  the  possibility  to  improve  beam  quality  in  presence  of 
strong  turbulence  -  with  D/r0  values  about  or  greater  then  10.  In  this  closed 
loop  system  we  achieved  1.5-6  increase  of  the  Strehl  parameter. 
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